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CnI ' Abstract 

O' 

t^^ ' The m-photon Jaynes-Cummings Hamiltonian is a natural generalization of 

o\ ■ 

O-i. the relevant operators for the time-dependent generalized m-photon Jaynes- 

^ ■ Cummings Hamiltonian. The dynamical equations for these operators are also 

^ ' given. These operators are needed and indeed are the basic building blocks for 



the much studied Jaynes-Cummings Hamiltonian. In this short note we give 



performing calculations in the context of the Maximum Entropy Formalism. 
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I. INTRODUCTION 

In a series of papers [|l]-^ , the generalized time-dependent Jaynes-Cummings Hamilto- 
nian in the context of Maximum Entropy Principle [MEP] and group theory based methods 
1^ was studied. In particular, in [Q] the MEP formalism was used to solve time- dependent 
N-level systems. A set of generalized Bloch equations, in terms of relevant operators was 
obtained and as an example the N = 2 case was solved. It was thus demonstrated in ||l| that 
the dynamics and thermodynamics of a two-level system coupled to a classical field can be 
fully described in the framework of MEP and group theory based methods. Further in [^ a 
time-dependent generalization of the JCM was studied and by showing that the initial con- 
ditions of the operators are determined by the MEP density matrix the authors were able 
to demonstrate that inclusion of temperature turns the problem into a thermodynamical 
one. An exact solution was also presented in the time independent case. Finally in Q] more 
detailed analysis of the three set of relevant operators was given. These set of operators 
are related to each other by isomorphisms which allowed the authors to consider the case of 
mixed initial conditions. 

As is well-known when we consider a quantum two-level system interacting with a single 
mode of quantized field one is led to the familiar Jaynes-Cummings Hamiltonian [JCH] [Q 
provided one is interested only in the difference of the population of the two levels. The JCH 
has been extensively used as a model Hamiltonian in fields such as quantum optics, nuclear 
magnetic resonance, and quantum electronics. A interesting study of the JCH is the periodic 
spontaneous collapse and revival due to quantum granularity of the field P] . In the rotating 
wave approximation [RWA], the JCH becomes solvable and it has been broadly used in the 
last years |f^-[T2 . 



The mean values of the field's population, correlation functions and nth-order coher- 



ence functions are of interest and useful in several applications. The MEP formalism allows 
us to describe a Hamiltonian system in terms of those, and only those, quantum operators 
relevant to the problem at hand. Thus, this formalism is suitable to study the Hamiltonian 
given in 0,^. In [^^ the population of each level and not their difference is considered 
therefore the resulting Hamiltonian is called a generalized time-dependent JCH. 

Our aim is to consider the m-photon generalized time- dependent JCH. The m-photon 
JCH describes the interaction of a two-level system with a single quantized mode of electro- 
magnetic field via m photon emission and absorption processes between the two-levels. A 
couple of forms of the m-photon JCH have been suggested |p!3| , p^ one of them being intensity 
dependent. This simply means that the coupling of this m-photon JCH is proportional to 
the square root of the number operator for the photons. As the validity of this effective 
Hamiltonian may be questionable under certain circumstances |J5[, we do not consider it 



in this paper. We consider the time- dependent generalized version of the simpler m-photon 
JCH in this paper. The question arises does at least one set of relevant operators exist 
for the m-photon JCM? The aims of this short note is to give such a set and to give the 
evolution equations for it. The layout of this paper is as follows. In the next section we 
recall some well-known results of the group theory based MEP formalism. In section two we 
give the relevant operators and the evolution equations for their expectation values i.e. we 
use generalized Ehrenfest theorem to obtain Bloch equations. 

II. SUMMARY OF THE MEP FORMALISM 

It is instructive to summarize the principal concepts of the MEP |T|-|3| Jl6| , [T7|| . Given 



the expectation values <0j > of the operators Oj, the statistical operator p{t) is defined by 

p(t) = exp[-Ao/-|:AAj, (1) 



where L is a natural number or infinity, and the L+1 Lagrange multiphers Xj, are determined 
to fulfill the set of constraints 

<0,>=Tr[p(t)6,], j = 0,l,...,L, (2) 

{Oq = J is the identity operator) and the normalization in order to maximize the entropy, 
defined (in units of the Boltzmann constant) by 

Sip) = -Tr[p\np]. (3) 

Eq. ID is a generalization of the more familiar density operator. For e.g. in open system, 
where we have Grand Canonical Ensemble there are two Lagrange multipliers, /? = ^y ^^^ 
fi are present, and we write the density operator as [|1^ 

p{t) = exp (pn{T, V, p)-pH + PpN) , (4) 

As is well-known the dynamics are governed by the time evolution of the statistical operator. 
The time evolution of the statistical operator is given by 

^^ft = [H{t),m]- (5) 

The essence of the MEP formalism in conjunction with the group theory method is to find 
the relevant operators entering Eq. 0) so as to guarantee not only that S is maximum, but 
also is a constant of motion. Introducing the natural logarithm of Eq. |I] into Eq. |^ it can be 
easily verified that the relevant operators are those that close a semi-Lie algebra under 
commutation with the Hamiltonian H, i.e. 

L 

[H{t),dj]=thY,g,,{t)d,. (6) 

Thus the relevant operators may be defined as those satisfying the above equation. Eq. ^ 
defines an L x L matrix G and constitutes the central requirement to be fulfilled by the 



operators entering in the density matrix. The Liouville Eq. |^ can be replaced by a set of 
coupled equations for the mean values of the relevant operators or the Lagrange multipliers 



as follows |T9[: 



^^|^ = -E^.,<0.>, j = 0, 1, ...,L, (7) 



4 = 



d\j J^ . 

-jr = l^^i9jh j = 0, 1,...,L. (8) 

^'' i=0 

In the MEP formalism, the mean value of the operators and the Lagrange multipliers belongs 
to dual spaces which are related by [|l^] 

<«'— W (^' 

III. THE RELEVANT OPERATORS AND EVOLUTION EQUATIONS 

The generalized time- dependent ni-photon JCH in the RWA takes the form 

H = Eib\hi + E2b\h2 + uoa^a + T{t) [-idHib\ + 7*626! a'f"^) , (10) 

(^ = 1), where 7 is the coupling constant between the system and the external field, Ej and 
u) are the energies of the levels and the field, respectively, a\ a, are boson operators, 6] and 
bj are fermion operators and T(t) is an arbitrary function of time. On setting m = 1 we 
recover the generalized time- dependent JCH oi Gruver et al., 0. 

Working in the context of the generalized time-dependent JCH Gruver et al., found 
that the relevant operators can be presented in three different but equivalent forms, each of 
them having different physical interpretations. These sets are connected via isomorphisms 
which allows one to go from one set into another. As noted in the advantage of this 
multiple representation comes from the fact that when partial information in any set is 



known, for instance the initial values only for some operators are known, it is possible to 
complete the missing information via the isomorphisms, if the complementary data in any 
other set is known, i.e. mixed initial conditions |T^. We find that working in the framework 
of the Hamiltonian llO] the same arguments go through. From the form of the Hamiltonian 



10| , we can guess by looking at the level's population and the structure of the interaction 



terms that a basic set of relevant operators satisfying Eq. ^ is 

N, = blh , (11) 

N2 = blh , (12) 

A = a^a , (13) 

I"^ = -fciHibl + Yhbla^"" , (14) 

F"" = i{-icrhih\ - 7*62^1at™) , (15) 

7V2,i = h\h2h{h^ . (16) 

These fundamental operators appear in the three possible sets of relevant operators outlined 
below. The above operators possess a simple physical interpretation. A'; is the number 
or population operator for level one, N2 is the number operator for level two. A is the 
familiar number operator for the photon/external field. The operator J™ may be considered 
as representing the interaction energy between the levels and the external the field. The 
particle's current between levels is governed by F™. A^2,i is the double occupation number 



operator. The operators pl] - |T^ and |TJ-|T5| can be considered as the m-photon quantum counter 
parts of the operators obtained for the semiclassical 2- level system studied in ref . |l| . 

The simplest set which closes a semi-Lie algebra with the Hamiltonian, Eq. ^ is found 
to be, 

iVf = {a^y Ni (a)", (17) 



Nl^ = (at)" iV2 (a)", (18) 

A" = (at)" A (a)", (19) 

jn,m ^ ^^t)n fm ^^y ^ ^20) 

^n,m _ ^^t^n ^m ^^^n ^ ^21) 

m^ = (at)" iV2,i (a)" , (22) 



n = 0, 1, . . .. For n = Eqs. p!7| - p2| reduce to the fundamental set of operators given in Eqs. |TT] 



through |T6|. This set of relevant operators is suitable for numerical simulation as it provides 
the simplest form of the system of differential equations for the evolution of their mean 
values. Eqs. p!7| - |2^ are the fundamental operators sandwiched between powers of creation at, 
and destruction a photon operators. This leads us to consider the operators with n > 1 as 
a measure of virtual transitions due to the absorption of more than one photon followed by 
the emission of the extra photons in a transition between the levels. 

It is important to emphasize, as pointed out in [0] that the assumption of RWA made 
at the beginning introduces in a natural way the set of correlation functions. For example if 
works out the commutator of the particle's current between levels, F"^, and the Hamiltonian 
one obtains correlation operators involving the population operators of levels with the field 
operators, see Eq. ^ below. Specializing Eq. ^ to the simplest case, n = 0,m = 1, we 
have averages of terms such as Nl = a'^Nia, which represent correlation between the field (a 
and at) and population of level one (^^"1). Clearly if one did not impose RWA, the structure 
of the algebra would be very different. Thus the algebra of the relevant operators carry 
the very blueprint of the physical assumptions we make. In short as is evident from Eq. ^ 
that the generator of the algebra, the Hamiltonian, generates a set of operators which are 
closely related to the physics of the problem and in this sense we mean that the set is 
physically relevant, thus it is not surprising that any physical assumptions that are put into 
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the Hainiltonian will be "transferred" to the operators which close the algebra with the 
Hamiltonian. 

The other two sets of relevant operators which satisfy Eq. are 



Oi {a)T{aT + {a)T{aTd. 



n=0 



and 



1 r 



Oi {a^aY + {a^aY Oi 



n=0 



Here Oi are the fundamental operators given by Eq. |TT] through Eq. ^. This once again 
generalizes the m = 1 case considered by 0. As pointed out in 0] the first set can be 
interpreted as the correlation functions between the fundamental operators and (a^)"(a)". 
{a))^{a)'^ are proportional to the nth-order coherence function of the field, 0. The operators 
included in the second set are proportional to the correlations between the fundamental 
operators and the energy of the field. 

The dynamical equations for the operators given in Eq. p!7| - p2| can be obtained using 
the Ehrenfest theorem (Eq. |^), and are given by 

d < N^ > - 

-^— = T{t) <F"''"> +T(t)['^Cim <F"-i'"> 

1X6 

'n—2,m - 



Jt 
Jt 



+"C2m(m-l) <F"-"'™> + 



.... +" C^m{m - l)(m - 2) <F"'™>], 

-T(t) <F"''">, 

-a </"'"*> +2|7|2T(t)[- <iVf+'"> 



(23) 
(24) 



(25) 



J ^ fn,m ^ 

'l^ ± = a<F^'"'>, (26) 

= +T{t) ["+^Cim < F"''" > +"+^C2m(m-l) <F"-i'"> +.... 

+"+^C„m(m-l)(m-2)... <F°''">], (27) 



0, (28) 



dt 



n = 0, 1, . . ., where a = E2 — E1 —u and "C^ = . _^-,,^, . Eqs. P3l - Pq are the exact dynamical 
evolution equations of the relevant operators for the generalized time- dependent m-photon 
JCH. They can be thought of as a kind of generalized Bloch equations for the quantum field 
case. As can be seen, the different order correlations are connected via the operators A^" 



and F"^'™ , see Eq. ^ and ^. As a check, if set m = 1 in Eqs. ^ through |28| we find that 
these reduce to the equations resulting from the single photon generalized time-dependent 
JCH considered by Gruver et al., 0, see their Eqs. 22-27. 

Comparing our fundamental set of operators, viz Eq. |11] through Eq. |1^ with the one 
given in we note that the operators which have not changed are Ai, A^i, Ai and A2,i. Thus 
one would expect that [although in our case we have more involved correlations between A^" 
and F"'*"], 

|< (a^)" Al (a)" > + < (at)" A2 (a)" > - < (at)"-iA(a)"-i >|°° , (29) 

and 

{<(at)«A2,i(a)">};;^^^, (30) 



are constants of the motion. Indeed we find this to be the case as can be seen from Eqs. 



28| . It is easy to see that [3^ holds by directly looking at Eq. |28|. To check the validity of ^9 

we first rewrite Eq. |2^ for n — 1 

(i < A"^^ > - 

= +r(t)pCim <F"-i'™> +"C2m(m - 1) <F"-2''™> +.... 

+"C„m(m-l)(m-2)... <F°'">] (31) 



Adding Eqs. ^ and ^ and subtracting ^ we immediately see that ^ holds. This provides 
an excellent check on our manipulations which are tedious. From the expressions ^ and ^ 
it follow that the particle's current between levels is equal to the photon's fiux. In the case, 
77, = we end-up with the conservation of the level's population. For n > we obtain a 
restriction for the correlations. 

It is clear from Eq. |2^ that the mean value of the operators are not be independent. 
This restricts the choice of the initial conditions. It is therefore necessary to choose a 
formalism which respects the restriction on the initial conditions. Following ref. [^ one 
may settle this issue by using the MEP density matrix given in Eq. [1|. As is known in the 
MEP formalism the mean values and the Lagrange multipliers live in the dual spaces. As 
mentioned in the Lagrange multipliers are numbers that can be freely chosen. Now once 
the restriction Eq. ^ on the mean values is implemented the equivalent restriction for the 
Lagrange multipliers are automatically satisfied when the density operator is diagonalized. 
In the MEP formalism used in [Q lack of knowledge on the mean value of one operator is 
equivalent to setting its Lagrange multiplier equal to zero. 

IV. CONCLUSIONS 

We have presented a generalized version of the m-photon JCH giving a description in 
terms of physical relevant operators. The temporal evolution equations have been worked 
out for one of the set which provides the simplest form of the system of differential equations 
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for the evolution of the mean values of the operators. Since an arbitrary function of time has 
been included, this formalism allows us to study the system even when the coupling is time 
dependent. Our work is a simple extension of the work of Gruver et al. to the m-photon 
case. One advantage, as mentioned in 0, of giving a description of the system in terms of 
the three sets of physical relevant operators is that it allows one to treat the case of mixed 
boundary conditions. 
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